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AN ANALYSIS OF ONCE-PER-REVOLUTION OSCILLATING AERODYNAMIC THRUST LOADS 
ON SINGLE-ROTATION PROPELLERS ON TRACTOR AIRPLANES AT ZERO YAW 1 

By Vernon L. Rogallo, Paul F. Yaggy, and John L. McCloud III 


SUMMARY 

A simplified procedure is shown jor calculating the once - 
per-r evolution oscillating aerodynamic thrust loads on pro- 
pellers of tractor airplanes at zero yaw. The only flow field 
information required for the application of the procedure is a 
knowledge of the upflow angles at the horizontal center line of 
the propeller disk. Methods are presented whereby these 
angles may be computed without recourse to experimental sur- 
vey of the flow field. 

The loads computed by the simplified procedure are com- 
pared unth those computed by a more rigorous method and the 
procedure is applied to several airplane configurations which 
are believed typical of current designs. The results are gener- 
ally satisfactory . 

INTRODUCTION 

An important consideration in the structural design of 
propellers is that which deals with vibratory stresses result- 
ing from the blade bending flatwise 1 cycle per propeller 
revolution (hereafter referred to as I P stresses). This 
flatwise bending of the propeller blade is due primarily to 
aerodynamic thrust loads oscillating 1 cycle per propeller 
revolution (hereafter referred to as 1 •P thrust loads) which 
are imposed on the blade as a result of propeller thrust- axis 
inclination and/or asymmetries of the flow fields in which 
the propellers operate. 

Methods have been developed whereby these 1 -P stresses 
have been computed satisfactorily for isolated propellers 
(e. g., refs. 1 and 2). These methods assumed a uniform 
flow field at the propeller disk, that is, the flow of air into 
the propeller disk was assumed to be at a uniform speed and 
stream angle at all points on the disk. However, for cases 
of propellers operating in nonuniform flow fields of wing- 
nacelle-fuselage combinations, the validity of these assump- 
tions was in question and it was not known to what degree 
the values thus computed would agree with those measured. 

In considering where errors might possibly arise in 
computing the 1-P stresses when the flow field is nonuni- 
form, it should be pointed out that prediction of the 1 P 
stresses requires a knowledge of the IT thrust loads and 
the structural properties of the propeller. In turn, methods 
of computing the 1 -P thrust loads require a knowledge of 
the flow-field characteristics, the blade physical properties 
(airfoil section, plan form, etc.) and the blade aerodynamic 
properties (two-dimensional). Finally, methods of com- 
puting the flow-field characteristics require taking proper 

1 Supersedes NACA TN 3395 entitled “On the Calculation of the VP Oscillating Aerodvn 
and Paul F. Yaggy, 1955. 


account of the influence of the various airplane components 
and their interference effects on each other. 

It was recognized that the blade structural, geometric, 
and aerodynamic characteristics did not depend upon the 
uniformity of the flow field; however, sources of error might 
possibly lie in the answers to one or more of the following 
questions. 

(1) If the IT aerodynamic thrust load is known, can the 
1 -P stresses be computed? 

(2) If (1) can be accomplished, can the IT thrust loads 
be computed if the flow-field characteristics are known? 

(3) If numbers (1) and (2) can be accomplished, can the 
flow field characteristics be computed for any given airplane 
configuration? 

To answer the above questions, an analysis of the methods 
employed to obtain 1 -P blade stresses, 1 -P thrust loads, and 
flow-field parameters was undertaken. An investigation 
of a propeller operating in the nonuniform flow field of a 
wing-nacelle-fuselage combination was made in the Ames 
40- by 80-foot wind tunnel. The data obtained from this 
investigation were reported in references 3 and 4. These 
data consisted of surveys of the flow field in the absence of 
the propeller, measurements of stresses on the propeller 
blades operating in this flow field, and surveys of the wake 
pressures behind the operating propeller. During the 
analysis of these data, a simplification was devised whereby 
the 1 -P thrust loads could be computed by a method re- 
quiring only a knowledge of the upflow angles at the hori- 
zontal center line of the propeller disk. Efforts were then 
directed toward developing methods whereby these upflow 
angles could be computed without recourse to experimental 
survey (these surveys had proved to be quite tedious and 
cumbersome). The methods which were developed were 
presented in references 5, 6, and 7. To evaluate the ade- 
quacy of these methods, it was necessary to obtain addi- 
tional experimental flow-field information on different 
wing-nacelle-fuselage combinations. The results of the 
surveys at the propeller planes of six model configurations 
tested in the Ames 40- by 80-foot wind tunnel were reported 
in reference 8. Mach number effects on the upflow angles 
were investigated to a Mach number of 0.92 by Messrs. 
Lopez and Dickson in the Ames 12-foot wind tunnel. These 
results were reported in reference 9. 

This report presents an analysis, based on the information 
in references 3 through 9, of the problem of computing 1 • P 
stresses for propellers operating in nonuniform flow fields. 

lie Loads on Single-Rotation Propellers in Pitch on Tractor Airplanes,” by Vernon L. Rogallo 
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NOTATION 

blade-section thrust coefficient, p ~ 2 jJi 

blade-section thrust coefficient due to blade-section 
section thrust due to lift 

lift force only, 
propeller diameter, ft 

up wash parameter, rate of change of upwash with 
angle of attack 

propeller rotational speed, revolutions per second 
except where otherwise noted 
propeller radius, ft 

distance along any radial line from the propeller 
thrust axis, ft 
blade-section thrust, lb 

free-stream velocity, feet per second except where 
otherwise noted 

component of the local velocity in the plane per- 
pendicular to a radial line (see fig. 1), ft/sec 
resultant velocity acting on any blade section, 
neglecting propeller-induced effects, ft/sec 
local velocity at any point (r,Q) at the propeller 
plane, ft/sec 

radial location of any blade section, 


1 P 


a 


exo 


fix 

8c 1 1 


Ac, 

4 *, 


V 

0 


e' 


the fundamental component (1 cycle per propeller 
revolution) of the oscillation of the blade-section 
thrust coefficient or of the vibratory blade stresses 
angle of attack of a wing or a body relative to the 
free-stream velocity measured at the plane of 
symmetry, deg 

geometric angle of attack of the thrust axis relative 
to the free-stream velocity, deg 
section blade angle at radial station x, deg 
amplitude of the \-P component of the variation oi 
incremental section thrust coefficient due to sec- 
tion lift force, further defined as the 1 P aerody- 
namic thrust load 

incremental section-thrust coefficient, c< Q ^(q = o 
incremental section-thrust coefficient due to section 
lift force, c t(Q — c , iQ=0 

angle of upwash, measured from a line parallel to 
the free-stream direction in a plane parallel to the 
model plane of symmetry, deg 
spanwise station from the longitudinal center line of 
a wing, semispans 

angle of outflow, measured from a line parallel to 
the thrust axis in a plane through the thrust axis 
(see fig. 1 ) 

angle at which the local velocity at any point on 



Fiernu 1 . —Geometric characteristics of the flow parameters 
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the propeller disk is inclined to the plane per- 
pendicular to the radial line through that point 
(see fig. 1 ) 

A sweep angle of the wing quarter-chord line, positive 
for sweepback 

A p effective wing sweep for compressible flow (A/?= 

tan --1 (tan A)/j8) 

p mass density of the air in the free stream, slugs/cu ft 

< p a phase angle between the position of the maximum 

l-P magnitude and the 12=90° position, deg 
i p angle of rotational flow (an apparent, not an actual 

rotation) measured from a line parallel to the 
thrust axis in a plane perpendicular to a radial 
line (see fig. 1), deg 

12 angular position about the thrust axis, measured 

counterclockwise from the upper vertical position 
as seen from the front, deg 

SUBSCRIPTS 


T total effects resulting from consideration of varia- 
tions in \p, 0, and acting simultaneously 
' 00 
v t 

0, y — effects resulting from consideration of variations in 

' CO 

y 

0 and -r—~ only 
* 00 

p effects resulting from consideration of variations in 

i/' only 


RESULTS AND DISCUSSION 

COMPUTATION OF THE IP STRESSES FROM KNOWN l-P AERODYNAMIC 

THRUST LOADS 

A method for computing the l-P vibratory stresses on a 
propeller blade using known l-P thrust loads has been de- 
veloped and presented in reference 10. Although this 
method has been proven to be adequate for isolated pro- 
pellers operating in uniform flow fields, its accuracy for 
propellers operating in nonuniform flow fields has not been 
established. An evaluation of this method was undertaken 
for a typical wing-nacelle-fuselage combination. Wind- 
tunnel tests of the combination were made in the Ames 
40- by 80-foot wind tunnel to provide the necessary ex- 
perimental data. The propeller stresses were measured 
on the blades by means of sti*ain gages. Simultaneously, 
the variations of the aerodynamic thrust loads were measured 
by means of total-head pressure rakes mounted in the pro- 
peller slip-stream as shown in figure 2. From the variations 
thus obtained, the l-P components of both the vibratory 
sti’esses and the oscillatory thrust loads were extracted by 
conventional wave-form analysis. 

A typical radial variation of the double amplitude of the 
l-P thrust load is shown for a blade angle of 20° at 0.75R 
in figure 3. By use of these measured l-P thrust loads, 
l-P stresses were calculated and are compared with the 
measured l-P stresses in figure 4. (Also shown in fig. 4 is 
a similar comparison for another blade angle.) Good agree- 
ment is indicated which is interpreted as verifying that the 



Figure 2. — Twin-engine airplane mounted in the Ames 40- by 80-foot 
wind tunnel showing propeller and wake survey instrumentation. 



Figure 3. — Measured radial variation of the maximum l-P incremental 
thrust coefficient; 1^ = 165 mph, n„=1250 rpm, a <y = 8°, p Q . 75 ~20°, 
four-blade propeller. 


l-P stresses in a propeller operating in a nonuniform flow 
field can be computed if the l-P thrust load is known. 

COMPUTATION OF THE I P AERODYNAMIC THRUST LOAD FROM KNOWN 
FLOW-FIELD CHARACTERISTICS 

A method utilizing strip analysis to compute the thrust 
load variation from known flow-field characteristics has been 
developed and is presented in Appendix A. An evaluation 
of the method was undertaken for propellers in nonuniform 
flow fields by attempting to compute the l-P thrust loads 
measured during the tests of the wing-nacelle-fuselage 
combination shown in figure 2. 

The flow-field characteristics required for this wing- 
nacelle-fuselage combination were obtained by surveying 
the flow field at the plane of the propeller in the absence of 
the propeller. A rake of eight directional pitot-static tubes 
(described in ref. 3) was employed for the survey. The 
flow field was defined by three parameters; the rotational 
flow angle, ip, the outflow angle, 0, and the velocity ratio, 
Vil^ooi all of which are defined in the Notation and shown 
in figure 1. Typical variations of these parameters with 
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Figure 4. — Comparison of computed and measured stress distributions 
for a four-blade propeller. 


angular position about the thrust axis are shown in the 
following sketch. 



From the measured flow-field parameters, the variation 
of the total aerodynamic thrust load was computed by the 
method described in Appendix A. In figure 5(a), the 
computed variation is compared with that measured by 
the thrust rake for a typical case at the 0.7 radius of the 
propeller. The adequacy of the method is demonstrated 
by the agreement of both the absolute magnitude and the 
wave forms of the variations. The phase angle between the 
variations is of no consequence to the problem of computing 
the 1 P stresses; it resulted primarily from rotation of the 
slipstream before it reached the thrust rake. 



(a) Variation of the total incremental thrust coefficient. 

(b) Variation of the I P and 2 P components. 

Figure 5. — Comparison of computed and measured section thrust- 
coefficient variations at the 0.7-radius station of the propeller; 
V m = 1 05 mph, a <?=8°, n v — 1250 rpm, p 0 . 7b = 20°, four-blade 
propeller. 

From the variations shown in figure 5(a), the 1-P and 
2 -P oscillatory components 2 were separated and are com- 
pared in figure 5(b). Again, aside from the phase shift, 
good agreement between measured and computed values is 
indicated. Similar accuracy was realized for all radial 
positions on the propeller blade as is indicated in figure 6 
by a comparison of the variation of the 1-P thrust load 
with radial position as measured and as computed. 

To obtain the accuracy indicated in figures 5 and 6 in 
computing the 1-P thrust loads, it was necessary to com- 
pute, using the measured flow-field data, the variation of the 
total aerodynamic thrust load at at least 1G angular posi- 
tions for each radial position. This, in turn, required an 
accurate knowledge of the flow-field parameters at these 
angular positions. If this information is available, and the 
limitations set forth in Appendix A are not exceeded, it is 
concluded that the 1-P thrust load can be computed from 
known flow-field characteristics. 


2 No other harmonic components of appreciable amplitude were found to be present. 
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Figure 6. — Comparison of the computed and measured distributions 
of the maximum variations of the 1*P incremental section thrust 
coefficient; 1^ = 165 mph, n p = 1250 rpm, a G ~ 8°, /3 0 . 75 = 20°, four- 
blade propeller. 

COMPUTATION OF THE FLOW-FIELD CHARACTERISTICS FOR ANY GIVEN 
AIRPLANE CONFIGURATION 

On the basis of the results given previously, it appears 
that predictions of the I P stresses depend on an accurate 
knowledge of certain flow-field characteristics which can be 
obtained by sufficiently extensive experimental surveys of 
the flow field at the propeller plane. These surveys are 
quite laborious and it would be desirable to be able to compute 
the required flow parameters. It was found that no known 
methods for computing all the flow-field parameters were 
available and attempts to develop such methods posed 
formidable problems whose solutions were not readily 
apparent. Therefore, rather than attempting to approach 
this problem directly, an attempt was made to consider the 
influence of each of the flow-field parameters on the 1 -P 
thrust load; the prime purpose being to discover simplifica- 
tions which could be accepted in determining the amplitude 
of the 1-P thrust load. 

The variation of the oscillating aerodynamic thrust load 
on a blade element was calculated by equation (A8), Ap- 
pendix A, using measured flow-field parameters at 16 angular 
positions. Three cases were investigated; namely, 

(1) The use of measured values of ^ with the assumption 
that 9= 0°, and Vi/V m = l.O 

(2) The use of measured values of 9 and V t /V „ with the 
assumption that ^=0° 

(3) The use of measured values of \p, 9, and Vj/V „ 

From the variations thus obtained, the 1-P components 

were extracted. 

From figure 1, it is seen that the first case shows the effect 
of variations in the direction of the vector V' and the second 
case shows the effect of variations in the magnitude of the 
vector V', while the third case shows the combined effects 
of both. Although the effects of these parameters cannot 
be considered to be wholly independent, it is believed that 
the first-order effects were clearly indicated since the summa- 
tion of the loads computed for cases (1) and (2) was nearly 
identical to those computed for case (3). 

The study included four configurations, 3 one nacelle and 

3 The configurations selected were believed to be representative of current design trends. 

A four-bladed, right-hand propeller (described in ref. 4) was assumed to be installed on each 
configuration 
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three wing-fuselage-nacelle combinations, which are shown 
in figure 7. The flow fields at the propeller plane were 
surveyed by use of a rake of directional pitot-static tubes. 

The 1-P components of the calculated oscillating aero- 
dynamic thrust load variations at the 0.7-radial station of 
the propeller, and the associated flow-field parameters are 
shown in figure 8. It can be seen that the variation of \p, 
the direction of V', produced a 1-P component, (1-P)*, 
which was always phased such that it peaked at angular 
positions of approximately 90° and 270°, while the variation 
of 9 and V jjV a, , the magnitude of V', also produced a sig- 
nificant 1-P component (1-P)#, vuv a the phase of which 
varied for different configurations. Further, since the total 
1-P component, (1 -P) T) is the summation of these two com- 
ponents as expressed above, the phase and amplitude of the 
total 1-P component curve are directly dependent upon the 
amplitude ratio and phase relation of the two components 
due to \f/, and to 9 and Vi/V m . For the configurations in- 
vestigated, (1-P)* was predominant in amplitude, its am- 
plitude being never less than 87 percent of the total 1-P 
amplitude. Although for some configurations (1 -P)e. vi/v 
had significant amplitude, it was phased such that its ad- 
dition to (1-P)* caused the amplitude of (1-P) r to be only 
slightly larger than that of (1-P)*. For other configurations, 
phasing of (1-P)*, viiva> was such that it could have in- 
creased the amplitude of (1 -P) T over that of (1-P)*. How- 
ever its amplitude relative to that of (1-P)* was too small 
to be of significance. Insofar as the configurations inves- 
tigated are typical of current design trends, it would seem 
that the amplitude of (1-P)* would be a good approxima- 
tion of the amplitude of (1-P) r , the value necessary for 
computing the 1-P stresses. It also is noted that (1-P)* 
always peaked at or ver} r near the 90° and 270° angular 
positions and that at these positions, 0 and VJV „ were 
generally 0° and 1.0, respectively. 

DEVELOPMENT OF A SIMPLIFIED PROCEDURE FOR DETERMINING THE 
AMPLITUDE OF THE I P AERODYNAMIC THRUST LOAD 

During the course of these studies, an important fact was 
discovered concerning the components of the oscillating 
aerodynamic load. For the cases investigated, no significant 
odd-order components above the fundamental, 1-P, were 
present. This fact makes possible the extraction of the 1-P 
component from the total variation by a simple operation 
which is independent of the number or amplitude of even- 
order components present. This simplicity stems from the 
fact that even-order harmonics repeat their magnitude and 
sense at half-cycle intervals of the fundamental, while odd- 
order components repeat their magnitude but are of opposite 
sense. If odd-order components above the fundamental 
are not present, one-half the arithmetic difference between 
values of the oscillating thrust load at any two angular 
positions 180° apart yields the exact magnitude of the 1-P 
component at those positions. Thus, computations at 
several angular positions around the propeller disk will 
define the exact variation of the 1-P component. 

The amplitude and phase relation of the 1-P component 
can be determined by applying the method described at any 
four equallj^ spaced angular positions. This is possible 
because (1) the variation of the oscillating aerodynamic 
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(a) Model A. 


Wing characteristics 




Wing characteristics 



(0 


(c) Model C. 



Wing characteristics 


Aspect ratio 

7.3 

Section normal to c/4 line 

64A010 

Area 

233ft 2 

Sweep of c/4 line 

40° 

Wing incidence 

3.4° 

Taper ratio 

0.49 

Scale 

0.41 


Propeller disk-y 



(d) Model D. 


Figure 7. — Geometric characteristics of the models. All dimensions in inches. 
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Jl-P>r 
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i i 



0 90 180 270 360 


Angular position, R, deg 


(b) Model B. 



Angular position, R, deg 


(c) Model C inboard nacelle. 


Figure 8. — The flow-field parameters and their effects on the 1-P variation of incremental section thrust coefficient for an isolated 

nacelle and several wing-fuselage-nacelle combinations. 
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Figure 8.— Concluded. 


load was periodic, (2) the components were all integer func- 
tions of one propeller revolution, and (3) no odd-order com- 
ponents above the fundamental were present. For example, 
at the specific angular positions 12= 0°, 90°, 180°, and 270°, 
expressions for determining the phase and amplitude of the 
1 P component are 

\ [C I( )n=90"-( c </)n=2-o»]= 5(, < ( cos ^ Ca) 

\ [( c <,)n=o°-( c ',)a-i80»]= 5c ‘i sin to (lb) 


where <p a is the phase angle between the position at which 
the maximum magnitude of the 1 -P component occurs and 
the 12=90° position. As noted previously, for the cases 
investigated, values of <p a were found to be quite small. 
Hence, it is believed plausible to make the approximation 
that 

8c 1 1 cos <p a ~8c tl (2) 

A simplified procedure for determining the amplitude of the 
1 P aerodynamic thrust load is now evident which requires 
only the solution of equation (la) with the foregoing approx- 
imation. 
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Figure 9. — Comparison of I P thrust loads as computed by the 
procedure of this report with those computed using the entire 
measured flow field for an isolated nacelle and several wing-fuselage- 
nacelle combinations; a?=0.7, I r ro = 165 mph, n D = 1250 rpm, 
/3 0 .7o=21.7°. 

The introduction of this simplified procedure reduces the 
flow field information required to the values of the param- 
eters at the 12=90° and 270° positions (the horizontal center 
line of the propeller disk). This fact and certain character- 
istics of the flow-field parameters preclude the need for 
experimental surveys of the flow field. It can be seen in 
figure 8 that, generally speaking, at the 12=90° and 270° 
positions 6 and VJV « are near the free-stream values of 0° 
and 1.0, respectively, as noted previously, while is at 
nearly its maximum and minimum values. Hence, 0 and 
Vi/Vo, can be approximated by the free-stream values. 
Values of i/' on the horizontal center line of the propeller disk 
may be computed by the methods which are presented and 
evaluated in Appendixes B, C, and D. Thus, the need for 
experimental surveys of the flow field has been obviated. 

To indicate the accuracy and general applicability of the 
procedure, 1 P thrust loads on a propeller computed by the 
procedure are compared in figure 9 with those presented 
in figure 8 which were computed by 16-point computations 
based on all measured flow parameters. As noted above, the 
values of 0° and 1.0 were assumed for 0 and VJV m when 
applying the simplified procedure. Comparisons in figure 
9 are at the 0.7-rad ial station of the propeller blade for all 
the test configurations. Comparisons of radial variations 
for two of the models are shown in figure 10 and, as further 
evaluation, the 1 -P blade stresses computed using these load 
variations are compared in figure 11. From these figures, 
it is seen that the simplified procedure yielded generally 
satisfactory results for these configurations. 4 

4 A convenient design factor, now in common use, for indicating the relative stress levels on 
propellers for various flight conditions is the 1-P load exciting parameter Aq (refs. 1 and 2). 
By definition, A is identical to the absolute values of ^ at the horizontal center line of the 
Propeller disk. 



Fraction of tip radius, x 


(a) Model B port propeller. 

(b) Model B starboard propeller. 

(c) Model C port inboard propeller. 

(d) Model C port outboard propeller. 


Figure 10. — Comparison of the radial variation of the I P thrust load 
computed by the procedure of this report with that computed using 
the entire measured flow field; U 00 = 16o mph, n p =1250 rpm, 
00.70=21.7°, 10°. 


IP vibratory stress, psi 
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(a) Model B. 


(b) Model C. 


Figure 1 1 .—Comparison of 1-P vibratory stresses resulting from I P thrust loads computed by the procedure of this report with jdiose 
resulting from I P thrust loads computed using the entire measured flow field; \ ^ = 165 mph, w p =1250 rpm, /Stwo— 21.7 , g— 10 . 


Evaluations of the methods for computing the values of 
\p at the horizontal center line of the propeller disk are 
presented in the appendixes. 

CONCLUDING REMARKS 

The simplified procedure presented herein enables the 
rapid calculation of the 1-P thrust loads on propellers of 
tractor airplanes in pitch at zero yaw, without the need for 
tedious experimental surveys of the flow field at the pro- 
peller plane. The simplifications of the procedure have 
reduced the necessary flow-field information to a minimum; 
namely, the upflow angles at the horizontal center line of the 
propeller disk. These angles may he obtained simply by 


the theoretical methods which are presented. Thereby, 
the need for experimental survey of the flow field is 
eliminated. 

The evaluations of the simplified procedure which are 
shown indicate that the 1-P thrust loads computed for the 
airplane configurations investigated are generally satis- 
factory. These configurations are believed to be typical of 
current designs and it is believed that equal accuracy could 
be expected for configurations generally similar to those 
investigated. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 21 , 1955 


APPENDIX A 

COMPUTATION OF BLADE-SECTION THRUST AND BLADE-SECTION NORMAL FORCE 


The 1 P stress problem is primarily clue to flatwise bending 
of the propeller blade resulting from 1 P oscillations in the 
blade normal force as the propeller turns through one 
revolution. The blade normal force at any angular and 
radial position is the sum of the components of the aero- 
dynamic loads (i. e., the thrust and torque forces) perpen- 
dicular to the chord line of the blade element. 

For the purposes of this investigation, it was desirable to 
express the normal force in as simple terms as possible. The 
normal force on a blade element, N b , may be expressed 
as 

N b =t cos Px+/q sin fi x (Al) 

v here 

t section thrust force 
j q section torque force 

Equation (Al) may be rewritten in coefficient form and 
both the section thrust and the section torque may be 
expressed in terms of the section thrust due to blade lift 
only, c tr 


c ^-c (i (l-^^)cos ft+c (i (l 


, blade-section drag 

where tan 7= : : — - 

blade-section lift 


tan 7\ 
cot p) 


tan (p-\- 7) 


sin ft. 
(A2) 


Equation (A2) can be rearranged as follows 




cos 


l f COS [(7 — cp) 


-2(ft— g?)]+cos (7 +<p) 


cos (7 — ^>)-f~cos (7 - j - v?) 


} 


(A3) 


When conventional propellers are operated at low for- 
ward velocities (i. e., conditions corresponding to take-off 
and early climb where high 1 -P excitations are encountered), 
the blade sections have high values of lift to drag ratio; 
hence 7 is quite small as is the angle of attack of the blade 
section which is equal to the quantity (ft— <p). It can be 
seen that if these terms are neglected in the above expres- 
sion, the term in the bracket becomes unity which is tanta- 
mount to saying that the drag force on a blade element con- 
tributes little to the blade element normal force. It is 
cautioned that such approximation may not be valid for 
all possible operating conditions. For simplicity, this 
assumption regarding the normal force was made for the 
work of this report and the normal force was expressed by 
the approximate relation 


CNb cos ft 


(A4) 


Steady-state theories, applied in strip analysis form, have 
been developed to compute the aerodynamic loads on pro- 
pellers operating in symmetrical flow fields where all the 
inflow is perpendicular to the propeller plane (e. g., see 


ref. 11). It appeared reasonable that application of t his 
type of analysis in stepwise computations might enable 
satisfactory computations of the values of c tl for propellers 
operating in distorted flow fields (thrust axis inclined or 
noninclined). The blade-element loads depend upon the 
angle of attack and velocity which the element experiences. 
The velocity involved is that in the plane tangent to the 
blade section path of rotation, V B in figure 12. For a rotat- 
ing blade element, the magnitude of V B is defined by the 
inflow velocity V' and the rotational velocity 1 m P Dx. If 
the interference effects of the blades on each other are 
ignored, the angle of attack of the blade element is the 
blade geometric setting ft minus the helix angle p 0 . 


Line parallel to thrust axis 



Figure 12 . — Velocity diagram of propeller blade section. 


The distorted flow-field case differs from the symmetrical 
case in that V' is not always perpendicular to the propeller 
plane nor of constant magnitude, but varies in direction and 
magnitude with angular position Q of the blade. In figure 
12, subscript 1 represents the symmetrical case while sub- 
scripts 2 and 3 represent the extremes of the variation for 
the distorted case. Consequently, it is apparent that these 
variations must be taken into account in the computations. 

The equation for the coefficient of thrust on a blade 
element fora propeller in a symmetrical flow field including 
blade interference effects is 


where 


C t = K7^X 3 - *** - 

57.3 


(cot ip — tan 7) 


(A5) 


ip ip o ~\~ OC 1 

p 0 tan" 1 (V m /Tr n p Dx) 

a t propeller-induced angle of inflow, deg 

k Goldstein correction factor for a finite number of 

blades 

. /blade-section drag\ 

7 tan 1 ( -ri — i : 7T7r~ ) 

\ blade-section hit J 
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For the distorted flow-field case, it is necessary to modify 
the equation to account for variations in both direction and 
magnitude of V' described above. Thus, 




an 7) / V s in. A 

a i V \ irn v Dx ) 


57 ' 3 ( COt 


(A6) 


for which it is important to note that 

V' cos p 


<p 0 =tan 


-i / V ' CQS ^ \ 

Vtt n p Dx—V' sin p/ 


In determining a* (the interference effects), the assump- 
tion is made that all the blades are operating at the same 
flow conditions as the blade being considered. 

The expression for c t may be simplified to consider only 
the section thrust due to blade lift, c tl , by ignoring the blade 
drag term, tan 7. Thus 




57.3 


cot (p / V' sin A 2 
< Y \ irn p Dx ) 


( c ° tH "5?y 


(A7) 


This is the expression used in the studies to calculate the 
values of c tl using all measured values of the flow-field par- 
ameters. When the values of 0° and 1.0 are substituted 
for 6 and V t /V « as in the simplified procedure of the report, 
V' becomes equal to V m . Hence, T r „ is substituted for V' 
in the last term of equation (A7) and in the expression 
for <p 0 . 


C [ . = K7T 3 

1 57 -< 


cot cp A _F 00 sin \p 
a i \ 2 \ 7 rtlpDx 


where 


(cot <P+~) 

, / V m cos p \ 

=tan 1 ( Tf — • — , ) 

\ t m p Dx—v co sin p) 


y (as) 


APPENDIX B 


METHODS FOR DETERMINING THE VALUES OF ^ AT THE HORIZONTAL CENTER LINE OF A PROPELLER DISK 


The value of ^ at the horizontal center line of a propeller 
disk is the summation of the geometric angle of attack of 
the thrust axis and the upwash angles induced by t lie 
various components of the airplane at the horizontal center 
line of the propeller disk, expressed as 

^ = OfgcometricT" € wingT" ^adjacent bodies T *bndy containing thrust axis 

(Bl) 

To determine these angles theoretically, the total angle of 
attack of each component must be known so that its con- 
tribution to can be computed. Each component total 
angle of attack, like is the result of a geometric angle of 
attack plus an induced angle resulting from the upwash 
of other components. It is clear that to obtain an exact 
total angle of attack for each component, an iterative 
process would be required. However, within engineering 
accuracy, it was found possible to avoid this process. 

In t lie case of the wing, it was found that could be 
adequately predicted by considering only an isolated wing 
at the geometric angle of attack; that is, the adjacent 
bodies influenced the wing total angle of attack so slightly 
that the resultant change in e w i„ g was. negligible. 

This conclusion is verified in figure 13 where e„i„ g com- 
puted by the method presented in Appendix C from measured 
span loading on a wing affected by a nacelle is compared 
to that computed for the isolated wing. The differences 
are seen to be negligible. 

In the case of determining the correct total angle of attack 
for a body, it was found that the effect of upwash from the 
wing and from other bodies generally could not be ignored. 
The wing upwash contribution to the angle of attack of a 
body can be computed by considering only the wing geo- 
metric angle of attack, since bodies have been shown to 
have little effect on the wing induced upwash. For a wing- 
body combination (i. c., a fuselage or single-engine airplane) 
the body effective angle of attack is taken as the arithmetic 
sum of the geometric angle of attack of the body and the 




(a) Spanwise distribution of loading coefficient. 

(b) Spanwise distribution of upwash. 

Figure 13. — Comparison of measured spanwise load distribution for 
an isolated wing and a wing-nacelle combination, and the resultant 
computed upwash distribution at a distance of one chord ahead of 
the wing leading edge; A/=0.18, Cl = 0.32. 
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value of wing upwash in the propeller plane at the axis of the 
body 

**efly = ^geometric “i~ € wing (B2) 

For the case of a configuration having more than one body 
(i. e., multiengine-type airplanes where the propulsive 
units are housed in nacelles outboard of a fuselage) the total 
angle of attack of a second body, for illustration assumed 
to be a nacelle, is expressed as 


comparisons are made in figure 15 for the wing-nacelle- 
fuselage combination described in reference 9 for a Mach 
number range from 0.25 to 0.92. In computing the upwash 
for this combination, account was taken of compressibility 
effects, except for the effects on the body-induced upwash, 
whose mechanics were not clearly understood. Good agree- 
ment between computed and measured upwash angles is 
indicated up to a Mach number of about 0.8. Compressi- 
bility effects are discussed more fully in Appendixes C and D. 


a eff n ^geometric 4“ € wing 4" € adjacent b: dy (B3) 

where the value of e wlng is the magnitude of wing upwash 
in the propeller plane at the nacelle axis. The value of 
^adjacent body at the nacelle axis is computed for the adjacent 
body at the effective angle of attack determined from equa- 
tion (B2). This ignores the fact that the adjacent body 
is under the influence of the nacelle; this second-order effect 
is believed to be negligible. 

While wing upwash is always important, the importance 
of considering upwash from adjacent bodies when obtaining 
a body total angle of attack varies widely depending on body 
proximity, relative size, and longitudinal position. For 
example, in the case of a multiengined airplane, the effect 
of one nacelle on the total angle of attack of the other can 
be shown to be small if the first nacelle lies behind the 
propeller plane of the second. Thus wing sweep would vary 
the importance of upwash from adjacent nacelles in finding 
the total angle of attack of a given nacelle. On the other 
hand, if the adjacent body were large compared to a nacelle, 
for example, a fuselage where the fuselage nose was down- 
stream of the propeller plane, its effect on the nacelle total 
angle of attack would be important. Generally speaking, 
any body which extends ahead of the propeller plane at 
which \p is being computed will have a significant effect on 
the total angle of attack of the body containing the thrust 
axis corresponding to that propeller plane; for example, 
usually fuselages and inboard nacelles both must be taken 
into account in finding the total angle of attack of the out- 
board nacelle in four-engined swept-wing configurations, 
whereas only the fuselage significantly affects the inboard 
nacelle. 

To evaluate the use of equation (Bl) and the methods 
which have been described, the variations of the upwash 
angle due to the various components and the curve repre- 
senting the summation of these and the geometric angle of 
attack of the thrust axis are shown in figure 14 for several 
wing-nacelle-fuselage combinations at a Mach number of 
0.22. Also shown is a comparison of these results with the 
measured values which indicates good agreement. Similar 


0.17c* 





(a) Model B. 

Figure 14. — Measured and predicted upflow angles and predicted 
upwash components at the horizontal center line of the propeller 
disks of several wing-fuselage-nacelle combinations; ao=10°, 
A/=0.22. 
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Figure 14. — Concluded. 
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Figure 15. — The effects of Mach number on the upwash at the horizontal center line of the propeller disks of a wing-fuselage-nacelle combination 

utilizing a swept wing. 


APPENDIX C 


A METHOD FOR COMPUTING THE WING-INDUCED UPWASH ANGLES AT THE HORIZONTAL CENTER LINE OF A PROPELLER 

DISK 


PREDICTION OF THE UPWASH ANGLE INDUCED IN THE EXTENDED WING- 
CHORD PLANE BY WINGS OF ARBITRARY PLAN FORM 

It was shown in reference 12 that for a known symmetrical 
load distribution on a wing, the wing-induced downwash 
could be computed by use of influence coefficients designated 
a vn coefficients which were presented therein. The same 
basic equations derived for the downwash computations are 
applicable for upwash computations, but new values of the 
a vn coefficients are required. 

Tlie upwash can be found at specific points in the extended 
chord plane between the wing tips and ahead of the wing 
leading edge. The general expression for the upwash is 
given as 

309174—57 3 


where 



(Cl) 


G n 


M 

b 

c 


dimensionless circulation ttt- identical to the load 

bV m 

CiC 

coefficient at span station n 
free-stream Mach number 

span of a wing measured perpendicular to the vertical 
plane of symmetry, ft 

local chord of a wing measured in a plane parallel to the 
vertical plane of symmetry, ft 
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Ci 

n 


section lift coefficient ofj a wing, 


section lift 

(ZcoC 


an integer defining a spanwise station on the wing 
quarter-chord line for which the value of circulation 
is defined 

free-stream dynamic pressure, lb/sq ft 
induced velocity, perpendicular to the mean chord line 
of the wing, positive for upwash, ft/sec 

longitudinal distance from the wing quarter-chord line, 
ft 

compressibility factor, VI 

an integer defining a specific point within the wing plan 
form for which the boundary condition of no flow 
through the wing is applied 

X 

t longitudinal coordinate semispans 

The a m coefficients are obtained from figure 10. The 
specific stations at which the upwash may be found corre- 
spond to r = l, 2, 3, 4 or jj„=cos (rx/8) =0.924, 0.707, 0.383 
and 0, respectively. For convenience, the preceding equa- 
tion is rewritten so that the upwash angle and loading distri- 
bution are in terms of unit lift coefficient. Tims 


w 


0 

V 


( w/Vjo \ 
\( 3 ('hlKar/ t 


'TO 




(C2) 


where 


$C L 4 irpb 2 /S 




f, lC - span loading coefficient 

t L,Cav 

S wing area, sq ft 

Can average chord of a wing, ft 


C L lift coefficient of a wing, 


wing lift 


K„ average of the ratios of the experimental section lift 

27T 

curve slope to the theoretical value all at the 
same Mach number 

Both a n and c,clC L c a v are affected by compressibility, since 
they have been determined as functions of the parameters 
P(b' 2 /S)/K„ v and A^. 

MODIFICATION OF THE METHOD TO PREDICT UPWASH ANGLE AT POSI- 
TIONS ABOVE AND BELOW THE EXTENDED WING-CHORD PLANE 

This modification is based on the assumption that the 
vertical variation of the wing-induced upwash angle is sim- 
ilar to that induced by a simple horseshoe vortex; that is, 
that a percentage ratio can be established between the two. 
For low Mach numbers, the bound portion of the horseshoe 
vortex has the same sweep as the wing quarter-chord line, 1 
and lies in the wing chord plane. For a value of jj and t// 3, 
the vertical variation of upwash angle due to a horseshoe 
vortex can be obtained from the following equation which 
corresponds to equation (34) of reference 13. 


Vl T /l 3 +(l~ tj) tan Asl 2 +(1 — nF+f 2 

r 2r > tan Ag—r/g “If (rffi tan fa l { \ r < 3+(1_,?) tau Vi ]} 

|_(2rj t an A g -T/<3) 2 +.C7cos 2 AjLV c os 2 A fl 7 J U+(l + ’l)J 1>L . 

V[t/( 3+(1 — n) tan A (J ] 2 +(1 + i?) 2 +V 2 

UXsd ta n Ag ~cA ;) ~r^ tan 2 A 8 -r/d) 2 +/ycos 2 aJ[(^ tan 


V0? tan A^— r/0) 2 +7? 2 +f 2 


(C3) 


where 

v lateral coordinate from the wing longitudinal center 
line, semispans 

f vertical coordinate from the wing chord plane, semi- 
spans 

From the variation thus obtained, the upwash angle at 
any position above or below the extended chord plane may 
be expressed as a percentage of the value at the chord plane. 


Based on the foregoing assumption, the value obtained on 
the extended wing-chord plane for a wing of arbitrary plan 
form may be multiplied by this percentage value to obtain 
the upwash angle induced by the wing at the same position 
above or below the chord plane. A typical example of the 
variations of the upwash angle expressed in percentage of 
the wing-chord plane value are shown in figure 17 for a 
horseshoe vortex swept back 40°. 

■ For high Mach numbers, the horseshoe vortex should have a sweep angle A/> given by 
tan Ae=(tan A)//S in accordance with the Prandtl-Glaucrt rule. 
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EFFECTS OF COMPRESSIBILITY ON THE WING-INDUCED UPWASH 

The Prandtl-Glauert rule, which accounts for the effects 
of compressibility, is directly applicable to the subject 
method. The basis for the corrections applied to the method 
as simply stated in the Prandtl-Glauert rule is that as the 
Mach number is increased, the span load distribution of a 
wing distorts as though its longitudinal dimensions were 
increasing as the ratio 1/0. Thus, the effects of Mach num- 
ber on a given wing can readily be considered by finding 
the span loading at zero Mach number of a properly distorted 
wing. 

Since the wing-induced upwash distribution is directly 
dependent upon the span loading, the effects of compressi- 
bility on this factor were investigated. A typical example of 
the compressibility effects on a span loading is shown in 
figure 18. The small differences which are shown have 
little significance in terms of the upwash induced in the 
extended wing-chord plane. Hence, the span load for zero 
Mach number can be used satisfactorily. 

It can be seen that increasing the longitudinal dimensions 
of the wing results in increasing the angle of sweep and 
increasing each local chord while leaving the span and taper 
unaffected. Thus, the distance from the quarter-chord line 
to a given point in the chord plane must also be increased 
by the ratio 1/0. The compressibility effects on the effective 
sweep angle are sizable. For example, increasing the Mach 
number from 0.33 to 0.90 changes the effective sweep angle 
of a 40° sweptback wing from 4lK° to 62%°. In terms of 
the wing-induced upwash, this effectively increases the dis- 
tance from a point on the extended wing-chord plane to the 
vortex core and thereby reduces the core's influence at that 
point. Hence, the values of r and A must be corrected to 
r/0 and A/* before entering the charts of the ci vn coefficients. 

It is obvious that for a given value of r/0, the compressi- 
bility effects are not constant for all sweep angles. To 
illustrate this point, typical examples of the upwash angle 
distribution, spanwise for constant values of r/0, are shown 
for wings with 0° and 40° of sweep. Figure 19 shows there 


is no variation of the upwash angle with Mach number for 
an unswept wing at a given r/0 since tan A is zero. However, 
for the wing swept back 40°, it is shown in figure 20 that an 
increase from a Mach number of 0.33 to 0.90 produced 
sizeable changes in the upwash angle distributions at a 
constant r/0. 

PROCEDURE FOR APPLYING THE METHOD TO OBTAIN THE WING- 

INDUCED UPWASH ANGLES AT THE HORIZONTAL CENTER LINE OF 

THE PROPELLER PLANE 

The first step is to obtain the span load distribution 
(e. g., by the methods of ref. 12); Mach number effects on 
the span load need not be considered as shown previously. 
Second, determine the r/0 variation and the value of A^ at 
the particular Mach number for the propeller plane locations 
in the extended wing-chord plane. The maximum and 
minimum values of r/0 will define the range over which 
computations must be made. Next, compute the values of 
c/Cl at the four control points on the span for constant values 
of r/0 over the range defined above. This is accomplished 
by substituting the values of the span load and the a vn 
coefficients derived from the charts in equation (C2). If 
necessary, these values should then be modified for dis- 
placements above and below the chord plane by means of 
equation (C3). The fairing of the resulting data points 
will yield the spanwise distributions of e/C L for constant 
values of r/0; a plot similar to figure 20. Using the values 
of r/0 previously determined, the propeller plane can be 
located on this plot as shown in figure 20 and the values 
of e/C L determined at any radial station. 

To illustrate the method and the results obtained, the 
variations of r/0 and elC L across the horizontal center line 
of the propeller disks of two representative airplane con- 
figurations are shown in figures 21 and 22 for Mach numbers 
of 0.33 and 0.90. The geometric characteristics of the com- 
ponents of these configurations were identical except that 
one wing panel was unswept and the other was swept 
back 40°. 


Influence coefficient, 
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Figure 16. — Influence coefficients, a vn , for symmetric span loading plotted as a function of the sweep parameter, A*, for various distances 

ahead of the wing quarter chord line, r/0. 
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Sweep parameter, A/ 3 , deg 



(d) v= 1, n = 4. 

(e) p=2, n= ] . 
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Figure 16. — Continued. 
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Figure 1G. — Continued 
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(h) v = 2, n= 4. 




Sweep parameter, Ag, deg 
(k) v = S, n=3. 


Figure 16. — Continued. 


Influence coefficient, 
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(m) y=4, n= 1 . 

(n) v = 4, n = 2. 

(o) v=4, n= 3. 


Figure 16. — Continued 
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Influence coefficient, 
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(p) */ = 4, n= 4. 
Figure 16. — Concluded. 



Figure 17. — Variation of upwash of a 40° swept-back horseshoe 
vortex with vertical position. 
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(a) A = 0°. 

(b) A = 40°. 

Figure 18. — A comparison of effects of compressibility on the spanwise 
loading coefficient for an unswept and a swept-back wing. 



Figure 19. — The spanwise distribution of upwash at constant values of 
t//? for an unswept wing at subsonic Mach numbers. 



Figure 20. — The effects of compressibility on the spanwise distribution 
of upwash at constant values of r/0 for a wing swept back 40°. 
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Figure 21. — The variation of the longitudinal parameter, r/0, at the 
horizontal center line of the propeller disks of two hypothetical 
airplanes for a low- and a high-speed flight condition. 


Figure 22. — The variation of the upwash angle at the horizontal 
center line of the propeller disks of two hypothetical airplanes for a 
low- and high-speed flight condition. 


APPENDIX D 

METHOD FOR COMPUTING THE BODY-INDUCED UPWASH ANGLES AT THE HORIZONTAL CENTER LINE OF A PROPELLER DISK 


COMPUTATION OF THE UPWASH ANGLES AT THE HORIZONTAL PLANE 
OF SYMMETRY OF A BODY 

A method was developed for calculating the upwash 
angles induced at the horizontal plane of symmetry of a 
finite body of revolution with radius which varies along the 
longitudinal axis. 

The induced upwash at the horizontal plane of symmetry 
of a body of revolution is a result primarily of the displace- 
ment of the transverse flow about the body. The transverse 
velocity W (see fig. 1(b)) normal to the body center line 
results from the bod} 7 ' being at an angle of attack to the 
free-stream velocity as shown in figure 1(b). The transverse 


flow about an infinitely long cylinder may be obtained by 
covering the longitudinal axis of the cylinder with doublets 
of moment per unit length, ju (ref. 14). The potential 
function r/<t> of a doublet element of strength n dx is 

d$=-r-^-, cos Q sin X d\ (Dl) 

47 rr 

where 

v f 

X tan -1 — 7 
x 

r f distance along any radial line from the body longitudinal 
axis, ft 
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x' longitudinal distance from any doublet element on the 
body longitudinal axis to a transverse plane containing 
the point at which the upwasli angle is to be computed, 
positive downstream from the doublet element, ft 

The potential function for a distribution of doublets 
along the longitudinal axis of a finite cylinder is obtained 
by integrating from the leading edge or most forward point 
(LE) to the trailing edge or most rearward point (TE) 
on the body. 


sin X d\ 


— cos 12 f x ™ 

471 -/'' Jx LK M 

Differentiating equation (1)2) with respect to 12 

sin 12 p™ . 

, /i sm X d\ 

4?r r JV R 


d if 


(D2) 


(D3) 


The expression for the velocity Vq tangential to a body of 
revolution and perpendicular to the transverse plane (see 
fig. 1(b)), is 


T y l (IQ 
Vn ~r' c/12 


(D4) 


Substituting the value of <7<t> from equation (D3) in (D4) 




sin 12 pTE 

4t^TJ* le 


n sin X d\ 


(Do) 


If a strength of h=2t R b 2 W is assigned to each doublet 
representing a cylinder of constant radius /?&, equation 
(D5) yields the tangential velocity. 

For a body of revolution whose radius varies with position 
along the longitudinal axis, it was considered that each 
transverse section was a section of an infinitely long cylinder 
with radius equal to that at the section. The doublet 
strength n was allowed to vary as 7?p For this case, the 
equation for the tangential velocity Vq becomes 


Vq 


sin 12 V m sin a U x te 


(D6) 


When a and e are small angles, the upwasli angle in the 
horizontal plane for 12=90° is (see fig. 1(b)) 


Fn a f* X TE 

«===*,% />Y sin X d\ (D7) 

v » A r ) Jx LE 

As has been shown by Munk and von Kerman, for calcula- 
tions of potential flow in the vicinity of the nose of a body, 
variations in the body radius need be considered only over 
the forward half-body. The body downstream of the 
midsection can be assumed to continue to infinity with a 
radius equal to that at the midsection. Therefore, the 
upper limit X TE can be assumed to be it. The integration 
is accomplished numerically to the midsection, after which 
the integral becomes a simple one because of the assumed 
constant radius. 


The expression in equation (D7) is for closed bodies ot 
revolution. For bodies through which there was air flow, 
in applying the equation the surface was considered to 
close across the openings as if there were no through flow. 

It is realized that the axial component of velocity is 
somewhat increased over the free-stream velocity because 
of the variation in body shape along the longitudinal axis. 
However, these increases are only a small percentage of 
V m and, therefore, they have been neglected. Similarly, 
a radial flow is induced by the body shape variation, but 
since it does not enter into the computation of e at the 
horizontal plane of symmetry, it has been ignored. 

Minor asymmetries in body shape downstream of the 
propeller plane were found to have little effect on the up- 
wash angles induced at the propeller plane. This will be 
demonstrated during the evaluation of the method. If the 
body was not symmetrical in the longitudinal vertical plane, 
it was found satisfactory to use the mean thickness line of 
the body when determining the effective angle of attack. 



Figure 23. — The variation of upwasli angle with vertical position in a 
transverse plane for an infinite cylinder. 
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EXTENSION OF THE METHOD TO POSITIONS ABOVE AND BELOW THE 
HORIZONTAL PLANE OF SYMMETRY 


For angular positions other than the horizontal plane of 
symmetry of a body, the up wash is no longer a function of 
the induced tangential velocity only, but is also affected by 
the induced radial velocity. To obtain the upwasli at posi- 
tions other than the horizontal plane of symmetry for a 
finite body of varying radius is a difficult task for which 
there is, at present, no adequate solution. On the other 
hand, the upwasli at any point about an infinite cylinder in 
incompressible flow may be obtained quite easily. The 
value of e induced by an infinite cylinder at any point de- 
fined by the lateral and vertical distances from the body 
center line, y b and z b , for small values of a may be deter- 
mined as 


(W^) MzW 1 
V m - \ [WR>)*+WRtf ) 2 J 


(D8) 


Variations of e per unit angle of attack versus lateral posi- 
tion y h for several vertical positions z b are shown in figure 23. 

EVALUATION OF THE METHOD 

The method was applied at the horizontal plane of sym- 
metry of three bodies where the horizontal center line of the 
propeller disk was coincident with the horizontal plane of 
symmetry. The models were tested in the Ames 40- by 80- 
foot wind tunnel at a Mach number of 0.18 to obtain experi- 
mental values of the up wash angles (ref. 5). The geometric 
characteristics of the models are shown in figure 24. Fig- 
ure 25 show r s comparisons of the experimental and computed 
upwasli angle distributions. Good agreement is apparent. 

A fourth body, which was not symmetrical in the longi- 
tudinal-vertical plane, was tested in the Ames 12-foot wind 
tunnel (ref. 9) through a Mach number range from 0.33 to 
0.92. The geometric characteristics of the body are shown 
in figure 26. Computations were made for this body assum- 
ing the plane of symmetry to be that of the forebody and 
ignoring the as 3 T mmetry downstream. The comparisons 
shown in figure 27 at the lower Mach numbers indicate that 
the asymmetry had little effect on the accuracy of the com- 
putations. The measured and computed values do not 
agree beyond a Mach number of 0.6. These differences are 
believed due to the onset of compressibility effects which 
have not been taken into account in the computations. 

An indirect evaluation of the extension of the method 
above and below the horizontal plane of symmetry is shown 
by the results of the computations of total upwasli angles 
for complete wing-nacelle-fuselage configurations as shown 
in Appendix B. 
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Figure 24. — Geometric characteristics of three nacelles tested at 
M = 0.18. All dimensions are in inches. 
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(a) Basic nacelle. 


(b) Nacelle with conical fairing. 


Figure 25. — Comparison of measured and computed variations of the angle of upwash at the horizontal center line of the propeller disk. 

M=0.18. 
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(c) Nacelle with conical spinner. 
Figure 25 . — Concluded. 
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27. — The effects of Mach number on the upwash 
be/da for an isolated nacelle. 
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Figure 26. — Geometric characteristics of nacelle tested from Mach 
number of 0.33 to 0.92. All dimensions are in inches. 
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